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In this paper we study an evolution of low-temperature thermodynamical quantities 
for an electron gas with a 5-function attraction as the system crosses over from weak- 
coupling (BCS-type) to strong-coupling (Bose-type) superconductivity in three and two 
dimensions. 



I. INTRODUCTION 

Recently there was remarkable progress in the study of the crossover from weak-coupling (BCS) super- 
conductivity to Bose-Einstein condensation (BEC) of tightly bound fermion pairs. Although this crossover 
was first addressed a long time ago the recent interest in this subject |l(| was initiated by sug- 

gestions of its relevance to High- Temperature Superconductors [jlT]]- p0| ). This crossover is a remarkable 
example of smooth evolution between two very different physical systems. A few of its features are: in 
the BCS superconductors there is one characteristic temperature T c where superconductivity disappears 
due to thermal decomposition of Cooper pairs. In contrast, the "Bose-type" superconductor possesses 
two characteristic temperatures, namely a critical temperature which in the extreme Bose regime does 
not depend on coupling strength and a substantially higher temperature of thermal pair decomposition. 
Also the BCS condensate is characterized by a coherence length which coincides with the size of a Cooper 
pair whereas in the Bose limit the system possesses two characteristic length scales, namely a coherence 
length which increases with increasing coupling strength and a Cooper pair size which obviously decreases 
with increasing coupling strength. 

In fact, this crossover is related to many other phenomena in various branches of physics. In particular 
it was generalized to the Chiral Gross-Neveu (GN) model pl]fj^] . This model is considered as a toy 
model for QCD. There are also numerous ongoing discussions of the possibility of the appearance of a 
pseudogap phase in the Nambu-Jona-Lasinio (NJL) model [^3-25 



In contrast to the particle physics, in the theory of superconductivity the BCS-BEC crossover has been 
studied in great detail. However, the main subject of the study was the evolution of critical temperature 
and non- Fermi-liquid properties above T c . Intuitive clarity of the physical picture of this crossover is 
obscured by the absence of simple analytical and asymtotic results for many other physical quantities 
(e.g. for the behavior of thermodynamic functions). The low-temperature asymptotics of thermodynamic 
functions of a BCS superconductor are well known and can be found in many textbooks. However to the 
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1 The Gross-Neveu and Nambu-Jona-Lasinio models are models of N-component fermions fields that serve for 
description of dynamic chiral symmetry breakdown in particle physics. Originally the NJL model was proposed 
in analogy to BCS theory. The NJL model may be treated by mean-field methods in the regime of large number 
of colors N c with 1/N C serving as a small parameter of the problem. However in a real world N c = 3 thus 
as it is not infinite, the system may possess [e.g. at finite temperature] a fluctuation-driven restoration of the 
chiral symmetry while preserving nonzero constituent quark mass. This would be a phenomenon similar to the 
symmetry restoration without pairbreaking in strong coupling or low carrier density superconductors. However 
in the case of the NJL model there are complications that do not allow this problem to be addressed rigorously 
by perturbative methods |p3]-p5[. Thus BCS-BEC crossover in superconductors may in some sense serve as a 
valuable toy model for QCD. 
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best of the author's knowledge there are no analogous expressions published for the thermodynamics of 
BCS-BEC crossover. This paper is an attempt to fill in this gap. 

In the model to be investigated in this paper, the crossover from BCS-type to Bose-type superconduc- 
tivity takes place either by varying the coupling strength, or by decreasing the carrier density. Asymtotic 
calculations in the BCS and BEC limits will be performed using the crossover parameter xq that we adopt 
from |)| . This parameter is directly related to the ratio of the chemical potential and gap function at zero 
temperature. It is a monotonic function of coupling strength and carrier density [to be seen in Figs. ^ 
and in two and three dimensions, respectively]. It is also a direct measure for the scattering length a s 
of the electron-electron interaction, at zero temperature in three dimensions [to be seen in Eq. (|TT|)], and 
of the binding energy of the electron-electron pairs in two dimensions [to be seen in Eq. ( |l3|) ]. 

The plan of the paper is following: we start by a brief outlining of the results of Refs. || , |To| ] for the gap 
A(0) and chemical potential at zero temperature in the entire crossover region. These are subsequently 
extended by equations for the low-temperature (T/A(T) << 1) behavior of the gap functions, chemical 
potential and of thermodynamic functions in two and three dimensions in the BCS and BEC limits up 
to the boundaries of the region — 1 < x$ < 1 . 

In appendix we also calculate the evolution of thermodynamic quantities near T c in weak and moder- 
ately strong coupling superconductors where mean-field methods are reliable. 



II. THE MODEL 



The Hamiltonian of our model is the typical BCS Hamiltonian in D dimensions (h = 1) 

# = E jd D xij\(x)^l(x)i; l (x)^(x), (1) 

where ?/><t(x) is the Fermi field operator, a —~\, | denotes the spin components, m is the effective fermionic 
mass, and g < the strength of an attractive potential g<5(x — x'). 

The mean-field equations for the gap parameter A and the chemical potential \x are obtained in the 
standard way from the equations: 

-I = lyJ-tanh^, (2) 
g V ^ 2E k 2T ' w 

n=-Y (l- & tanh^V (3) 

k v ' 

where the sum runs over all wave vectors k, N is the total number of fermions, V the volume of the 
system, and E^ = + A 2 with £k ~ k 2 /2m — are the energies of single-particle excitations. 

Changing the sum over k to an integral over £ and over the directions of k, on which the integrand 
does not depend, we arrive in three dimensions at the gap equation: 



l^sfa; tanh^^, (4) 

9 3 J- It 2T W 

where the constant K3 = m ' /V2it 2 has dimension energy ' /volume. In two-dimensions, the density 
of states is constant, and the gap equation becomes 



— = K2 [ d£ ta.nVi — y "j" (5) 

with a constant K2 = m/27r of dimension energy -1 /two- volume. In two dimensions, the particle number 
in Eq. (||) can be integrated with the result: 
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2tt 



vV 2 + A 2 + ( u + 2Tlog 



1 + exp 



sJn 2 + A 2 



T 



(G) 



The (5-function potential produces an artificial divergence and requires regularization. A BCS super- 
conductor possesses a natural cutoff supplied by the Debye frequency ujjy. For the crossover problem to be 
treated here this is no longer a useful quantity, since in the strong-coupling limit all fermions participate 
in the interaction, not only those in a thin shell of width u>d around the Fermi surface. To be applicable 
in this regime, we renormalize the gap equation in three dimensions with the help of the experimentally 
observable s-wave scattering length a s , for which the low-energy limit of the two-body scattering process 
gives an equally divergent expression: 



1 1 \ m 
g + V ^k 2 



Eliminating g from (Q) and (||) we obtain a renormalized gap equation 



Air a* 



V ^ 



1 E k m 

tann tt 

2E k 2T k 2 



(7) 



(8) 



in which l/fcjra s plays the role of a dimensionless coupling constant which monotonically increases from 
— oo to oo as the bare coupling constant g runs from small (BCS limit) to large values (BE limit). This 
equation is to be solved simultaneously with (^|). In the BCS limit, the chemical potential [i does 
not differ much from the Fermi energy e^, whereas with increasing interaction strength, the distribution 
function broadens and /i decreases, and in the BE limit, on the other hand we have tightly bound 
pairs and nondcgcnerate fermions with a large negative chemical potential 3> T. 

In three dimensions at T = 0, equations (|^), (J^j were solved analytically in entire crossover region in 
|l to obtain A and /i as functions of crossover parameter l/kpa s . The results are 



A 

JL 



1 



[xoh(xo) + Hxo)} 2 ^ 3 

jJL A Xq 

A e F 
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Xoh(x ) - Ii(x a ) 



kFCl s 



irlxohixol + hixo)] 1 / 3 ' 



(9) 
(10) 

(11) 



with the functions 



n OO 

h(x ) = / dx 
Jo 



{x A - 2x x 2 +x 2 + l) 3 / : 



= (l + *g)V*J5(- «)- 



Ax 2 (l 



.3.2)1/4 F ( 2 '«). 



1 f°° 

h(x ) = - J dx 



1 1 7T 

(a; 4 - 2x x 2 + x 2 + l) 1 / 2 ~ 2(1 + x 2 ) 1 ^ V^' 



2 _ x\ 2 _ k 2 1 _ fl 2 _ y/1 + X 2 , + X 

K ~ (1 + ^)1/2^ " ' Xo ~A' Xl ~ 2 ' 

where E(j,k) and F(j,k) are the usual elliptic integrals. The quantities (||) and ( |l(i| ) are plotted as 
functions of the crossover parameter xq = /i/A in Fig. [l| 
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FIG. 1. Gap function A and chemical potential fi at zero temperature as functions of xo = A*/A in three 
dimensions. 

In two dimensions, a nonzero bound state energy eo exists for any coupling strength. The cutoff can 
therefore be eliminated by subtracting from the two-dimensional zero-temperature gap equation 

1 1 1 m f°° , 1 

-- = — > — = — / 12 

where z — k 2 /2mA — xq, the bound-state equation 

g V^k 2 /m + e 2tt J_ Xo 2z + e Q /A + 2x 

After performing the elementary integrals, one finds : e / A = y/\ + x 2 ] — xq. From Eq. (§) we see that 
at zero temperature, gap and chemical potential are related to xq by JlO|,p| 

- = — 7m. <"> 

JL = 2*o_ 

The two relations are plotted in Fig. ^. From ( |l4"| ) we find the dependence of the ratio co/ef on the 
crossover parameter xq: 



£0 = 2 v^+J^o (16) 

£F y/1 +X 2 + X 
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FIG. 2. Gap function A and chemical potential fi at zero temperature as functions of xq in two dimensions. 



III. THE CROSSOVER FROM BCS TO STRONG-COUPLING SUPERCONDUCTIVITY 

NEAR ZERO TEMPERATURE 

A. The range of validity of mean-field approximation 



We have extended the above relations to non-zero temperature in the BCS and BEC limits. We use 
a mean-field approximation which, in the range of its validity, allows a nice simple description of the 
crossover. It is known that in BCS limit the mean-field approximation is valid for the description of 
the entire superconductive region. In fact the system crosses over from BCS to Bose condensate in a 
very narrow region —1 < xo < 1 (see e.g. [0) so the BCS approximation is adequate for description 
of the entire superconductive region in the range Xo > 1- This is in contrast to BEC limit where the 
superconductive transition is governed by the thermal excitation of collective modes rather than thermal 
decomposition of individual particles, the effect which is missed in the mean-filed approximation. So in 
the BEC limit we must be restricted in the following discussion to the temperatures much lower than the 
temperature of superconductive transition: < T « Tbosc, where Tbosc = [n/2£(3/2)] 2 / 3 7r/m [see e.g. 

§|]. 

A special care should be taken when a mean-field approximation is used for the description of a two- 
dimensional system. It is well-known, that a straightforward calculations of the next-to-leading order 
corrections gives that the gap should be exactly zero at any finite temperature in 2+1 dimensions. This 
is due to the fact that, according to standard dimensional reduction arguments, a 2+1 -dimensional 
system is effectively two dimensional at high temperature and in two dimensions there is no Goldstone 
boson as it is articulated in Coleman-Mermin-Wagner-Hohenberg theorem p(| . However as it was shown 
by Witten [^7j (see also |2^,^,|l^,^|,^2) ) a direct calculation of corrections misses essential physics of a 
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two-dimensional problem and a mean-field results still have sense under certain conditions in 2D. The 
situation is following: the fluctuations can be made arbitrarily weak by decreasing temperature in 2+1 
dimensions (or e.g. increasing N in 2D zero temperature calculations in p7j , ^l| , p4| ) and then employing 
a "modulus-phase variables" A = |A|e l61 one can observe that the system possesses a very well-defined 
"mexican hat" effective potential that determines the gap modulus. Due to phase fluctuations in the 
degenerate valley of the potential, the average of the complex gap function is zero, however there exist 
a gap locally (i.e. in some sense the system in its low energy domain degenerates to a nonlinear sigma 
model). In 2+1 dimensions, as the temperature approaches zero the thermal fluctuations in the degenerate 
valley of the effective potential gradually vanish and at T — a local gap becomes a real gap. Existence 
of the local gap modulus A does not contradict the Coleman-Mermin-Wagner-Hohcnberg theorem since 
A is neutral under U(l) transformations. Thus at low temperature in 2+1 dimensions there appears an 
"almost" Goldstone boson that becomes a real Goldstone boson at exactly zero temperature. At low 
temperatures we can neglect contribution from weak phase fluctuation at temperatures much lower than 
the temperature of the Kosterlitz-Thouless transition, which in the weak coupling limit approximately 
coincides with the temperature of pair formation and in strong-coupling limit reaches a plateau value: 
Tkt — £f/8 = irn/8m, where = (ir/Tn)n is the Fermi energy of free fermions with the carrier density 
n and mass m . 

Let us now outline how finite temperature quantities will be derived. We start with inspection of the 
gap A behavior at low temperature. A(T) can be expressed with the help of the analytic zero-temperature 
results from For example in three dimension we can write an equation that relates the gap at zero 
temperature with the gap at finite temperature at a fixed crossover parameter xq (which enters implicitly 
A(0)): 

'dg ^ = fV ^ tanh^^. (17) 
+ A(0) 2 ./_„ + A(T) 2 2T 

This equation should be solved self-consistently with eq. (||) . We must observe that the above system of 
equations, in contrast to zero-temperature calculations has the temperature as an extra parameter. Thus 
we need to specify in what model we are studying temperature effects with two most general options being: 
(i) "fixed chemical potential model" - i.e. when we do not fix the carrier density but assume the presence 
of a reservoir which provides us with a temperature- independent chemical potential /i = /u(l/fcj?a s ; T = 0) 
or (ii) "fixed carrier" density model - i.e. when n is fixed and (i varies with temperature. We choose 
the first option since such a fixed-/x model is the most convenient for deriving asymptotic results for the 
finite-temperature behavior of the system. However the choice of the model is not principal since from 
(|J) it is seen that in a "fixed density" regime the chemical potential in BCS and BEC limits depends 
very weakly on temperature in comparison with the dependence on the coupling strength in the regimes 
where our formulas are valid. Moreover the temperature dependence of the chemical potential in the 
entire crossover region was calculated numerically in 2D in Jlif] and it is small comparing to its variation 
with coupling strength in the temperature range < T < Tkt- So the results are similar in both cases 
of the presence of a particle reservoir that provides temperature independent chemical potential when 
carrier density varies or in the model of temperature independent carrier density. 

In our calculation we use xo as the most convenient crossover parameter, since it depends via the simple 
relation xq = /i/A on the chemical potential which can be measured rather directly experimentally 
The parameter xq ranges from — oo in the strong-coupling (Bose-Einstein) limit to oo in the weak- 
coupling (BCS) limit. The relation between xo and the inverse reduced coupling strength between the 
electrons l/fci?a s is plotted for three-dimensional system in Fig. [}| The corresponding relation (p~6| ) in 
two dimensions between xq and the bound state energy eo of the electron pairs is plotted on Fig. 
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B. The gap functions 



Let us now turn to the region near zero temperature, where the mean-field approximation gives an 
exact results for the gap function. From ( |l7| ) we extract the asymptotic behavior in the three-dimensional 
case at low T having T/A(T) as a small parameter. For this purpose let us first rewrite ( fl7l ) in three 
dimensions as: 



+ a(o) 2 y/e+w? 



1 — tanh 



ye + a(t)^ 

2T 



(18) 



In the region xo > 1 one can assume density of states to be roughly constant, since the integrand of 
is peaked in the narrow region near f = 0. Then the integral in the r.h.s of (18) can be writen as 

ve 2 + A(T)2" 



/OO 
- fJ, 



Ve+wr- 



: exp 



T 



(19) 



This integral can be calculated e.g. by changing the integration varaible to y = \J (£/A) 2 + 1 — 1 and 
expanding the factor before exponent. The result is expressed via the error function and we arrive at the 
following expression for the behavior of the gap function near T = in the region xq > 1 



A(T) = A(0)-A(0)^/- 1 



T 



A(0) 



exp 



A(0) 
T 



1 + erf 



y /x* + l- l 
T/A(0) 



(20) 



where erf (x) is the error function. Since the density of states is nearly constant in this limit, the same 
equation holds in two-dimensions — apart from a modified gap A(0) given by (|l4|). 

In the weak-coupling limit, xq = /i/A(0) tends to infinity, and the expression above approaches expo- 
nentially fast the well-known BCS result: 



A(T) = A(0) - [2 7 rA(0)T] 1 / 2 exp 



A(0) 
T 



(21) 



For strong couplings with xq < — 1 the integration range, in contrast to the above case, does not include 
the point £ = 0, so we write the r.h.s of (pit) as: 



2VA(0) 



v/e/A(0)+x o 



-« y/wmr + (wf 



: exp 



y/(g/A(0)) 2 + (A(T)/A(0)) 2 
T/A(0) 



VA(0) 



(22) 



This integral can be calculated e.g. by changing the integration variable to v = J (^/A(0)) 2 + 1 — \J x^ + 1 
and expanding the factor before exponent so that the above integral is reduced to: 



2 V /A(0) 



Fo 



exp 





/>cx: 




T/A(0) 


/ exp 
Jo 


T 



Thus, we find the behavior of the gap function at xq < — 1: 



A(T) = A(0) - 



T 3/2 



^V^o a/A(0) 



exp 



\/vdv 



yV + A 2 (0) 
T 



(23) 



(24) 



From Eq. (|2J) we see that near T = the gap A(T) tends in the strong-coupling limit exponentially to 
A(0), forming plateau near T = 0. 
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In two dimensions we arrive at similar result: an exponentially growing plateau near T = in the 
strong coupling limit: 



A(T)=m .m El (v^±zy (25) 

where Ei is the exponential integral Ei{z) = J z °° e~ t /t dt. The plateau formation in the low temperature 
solution for the gap function in a strong-coupling regime (see Fig. 5) is an indication of the suppression 
of a pairbreaking effects in low temperature domain. One can however observe that in contrast to BCS 
regime in the BEC regime the suppression comes from a large negative chemical potential in the exponent 
H»|A|. 




FIG. 5. "Plateau" formation of the gap function at low temperature in three dimensions. In the BCS limit 
the thermal fluctuations start breaking fermion pairs already at low temperatures which results in the absence 
of plateau of A(T)/A(0). In the BEC limit in contrast we have a plateau near zero temperature. Inclusion 
of Gaussian corrections in the BEC limit leads to depletion of the average order parameter due to deple- 
tion of macroscopic occupation of zero momentum level. The order parameter in this limit becomes zero at 
T c = [n/2C(3/2)] 2/3 7r/m « 0.218eF. This region is not plotted on the above figure. 



C. The thermodynamics 

Let us now calculate thermodynamical quantities near T = 0. For the thermodynamic Gibbs potential 
n(T, pi, V) we have 



A 2 J&, + A 2 

tanh v ^„ 2Tlog 



2T 



2 cosh A^±^ 
2T 



+ 6 



(26) 



Here and in the sequel in this section, A(0) will be replaced by A. In three dimensions, Eq. (|2g) turns 
into the 
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n 



«3 / d£y/£ + [I 



: tanh 



2T 



2Tlog 2 cosh 



2T 



In two dimensions, we obtain instead: 



V 



= k 2 / d£ 



A 2 



2^ 



: tanh 



2T 



2Tlog 2 cosh 



2T 



(27) 



(28) 



We regularize the thermodynamic potential fl s of the condensate subtracting O n = fi(A = 0). At T = 
and for weak couplings this is found to depend on \x and A as follows: 



v 7 = 



V 



A 2 1 



1 



( , -n\n\ - -^vV + A 2 



In the BCS limit (xo — > oo) this reduces to the well-known result 



— = k 3 Vm 



A 2 
2 



In two dimensions, we have a formula valid for any strength of coupling or carrier density: 



V 



= K 2 



A 2 1 



H- W^ 2 + A 2 



In the opposite limit of strong couplings, we find in three dimensions: 

£=-^ 3 A^(-*o)- 3 / 2 . 



(29) 



(30) 



(31) 



(32) 



The gap A(0) has by Eq. (11) the strong-coupling limit A(0) ~ e_F[16/37r] 3 / 2 |a;o| 1 ^ 3 yielding the large-xo 
behavior 



5/2 7T /16\ 



15/4 



-2/3 



(33) 



In two dimensions, we substitute the gap function A of Eq. ([u]), into the thermodynamic potential (|3l|), 
and obtain for strong couplings an interesting result: 



7- 

Let us now turn to the entropy. In three dimensions near T = it is given for weak couplings by: 
S 



V 



2ttA 3 



T 



■ exp 



1 + erf 



V4 + i- i 

T/A 



2/iexp 



yV + A 2 

T 



(34) 



(35) 



For ///A — > oo, this reduces correctly to the BCS result: 



S r 

y = *3VM 



8ttA 3 



T 



ex P I - y 



(36) 



In two dimensions, the result is similar apart from the fact that A = A(0) should be given by Eq. ( |l4| ) 
and Ksy/Ji should be replaced by K2- In the strong-coupling limit where ///A <C — 1, we have for the 
entropy in three dimensions 
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s_ 

V 



/ 2xA2 ( V^ 2 + & 2 \ 

K 3— T ' + A 2 exp I I 



(37) 



and in two dimensions: 



S _ / 

— = -2^ exp I 



(38) 



From the entropy, we easily derive the heat capacity at a constant volume cy . In three dimensions it is 
given near T = for weak and moderate couplings by 



A 



cv = n 3y /Jl\/2irA 3 { exp 



1 + erf 



T/A 



(39) 



reducing in the limit xq — ► oo to the BCS result 



2A 



CV = K 3 y/JlV 27rA3 ^57J ex P 



(40) 



In two dimensions, the weak-coupling behavior is the same with the factor k^^/Jl replaced by k,i while 
the strong-coupling behavior in three dimensions is 



cy 



= *3^T-VV + A *)exp(-V?^ 



(41) 



and in two dimensions 



M 2 / 
cy = 2k 2 — exp I 



(42) 



IV. CONCLUSION 



The BCS-BEC crossover and related phenomena is now a subject of increasing interest in many branches 
of physics. In this paper we have presented simple asymtotical results for one of the important aspects of 
this phenomenon in superconductors. Namely we studied the behavior of thermodynamic functions when 
a system crosses over from the BCS superconductivity to the Bose condensate of tightly bound pairs. 
The first observation that could be made analyzing the thermodynamic expressions is that the crossover 
takes place in a quite a narrow region 1 > xq > — 1 which is in agreement with previous studies. 

While the crossover is known to be continuous, the thermodynamic functions assume an exotic form 
in the regime of strong attraction or low carrier density. The interesting feature is that a modified gap 
function \J A 2 + y? enters all the expressions in the BEC limit in a similar way as the gap function A 
enters thermodynamic expressions in the BCS limit. Since in the BEC regime A << — fi, the chemical 
potential is in some sense playing the same role as the gap function in the BCS theory. This result, 
that may seem surprising from the point of view of the BCS theory, has the following roots: The gap in 
the spectrum of single-particle excitations has a special feature 0], 0, @| when the chemical potential 
changes its sign. The sign change occurs at the minimum of the Bogoliubov quasiparticle energy 
where this energy defines the gap energy in the quasiparticle spectrum: 

£ gap = min(^ + A 2 ) 1/2 . (43) 



11 



Thus, for positive chemical potential, the gap energy is given directly by the gap function A, whereas for 
negative chemical potential, it is larger than that: 



(M 2 



A 

- A 2 ) 1 / 2 



for 
for 



(i > 0, 
\i < 0. 



(44) 



That means that for the positive values of the chemical potential the gap function has the same meaning 
as the standard superconductive gap. In contrast, in the Bose limit, when chemical potential goes below 
the bottom of the band, one should distinguish between the two notions of the "gap". At first, there is 
the order parameter A that reaches zero at T c [ see e.g. ||]. Also, in this limit, there exist the modified 
gap function ( [III ) that tends to \fi\ = E^/2, where E\, — l/ma 2 s thus being proportional to the binding 
energy of the pairs. So when one is concerned about e.g. low- temperature thermodynamic functions one 
should understand under the "gap" basically the chemical potential rather than the order parameter A. 
Other well-known and to some extend similar features arising when a system crosses over from BCS to 
BEC limit are: (i) the appearance of two characteristic temperatures of different physical origin T c and 
T* in BEC limit g,| in contrast to only one characteristic temperature T c in the BCS limit; (ii) the 
separation of notions of the Cooper pair size and coherence length in strong coupling condensate which 
in contrast coincide in the BCS limit j|,|]j7]]. 

The above results for thermodynamical quantities once again illustrate the interesting circumstance of 
the essentially richer nature of a strong coupling condensate comparing to a BCS condensate that helps 
to understand better the phenomenon of symmetry breakdown in Fermi systems. 
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APPENDIX A: THERMODYNAMICS NEAR CRITICAL TEMPERATURE OF 
SUPERCONDUCTORS WITH MODERATELY STRONG COUPLING STRENGTH 

In this Appendix we present expressions for the behavior of the gap and the thermodynamic functions 
near T c in superconductors with moderately strong coupling strength. We employ a mean-field approx- 
imation and thus our expressions are not valid in the crossover region where one should incorporate 
fluctuation corrections in the regime near T c . 

In the regime of weak or moderately strong attraction and also for moderately low carrier density the 
gap function behaves near T c as follows: 



A(T) 

2T r . 



In the limit /x/2T c 



T 

1 



1 + tanh 

2T r 



1 



1 , /I 

- — ; -Tr tanh 

( M /2T C ) 2 2T c 



_fi/2T ( 

oo this tends to the BCS result 



2 u 
1 H — arctan 

7T TlTr 



(Al) 



(A.2) 



The ratio of zero temperature gap to the critical temperature has the following dependence on crossover 
parameter close to the BCS regime: 
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A(0) _ tt / A(0) 



-1/2 



= - 1- 



1 

4^ 



-1/2 



(A.3) 



Whereas the critical temperature: has the following dependence on xq\ 

T r e 1 ( 1 3 



€F 71" V^O oXq 

The thermodynamic potential in three dimensions in the weak-coupling regime near T c is given by 



(A.4) 



(T c - T)A 2 



+ 



A 4 1 



2T, 



1 + tanh — 
2T, 



+ 



(A.5) 



4Vm/2T c ( M /2T c )2 tanh 2T C 



4 (2T C ) 2 

In the BCS limit, this reduces to the well-known formula: 



+ 



2 n 
1 H — arctan — — - 

7T TT1, 



^ = -« 3 VmA 2 ( i 



T 



1 A 2 \ 



t c 2tt 2 t 2 y 



The entropy behaves near T c in three dimensions in the weak-coupling regime like 

1 + tanh 



S s S S n — A 2 

v = = 



2T, 



with the BCS limit 



S s ^A 2 



(A.6) 



(A.7) 



(A.8) 



In order to derive the specific heat we must take into account the temperature dependence of the gap. 
In three dimensions, we find in the weak-coupling region near T c : 



a 



1 + tanh — 
2T r - 



1 1 . , I 1 

— ; ; — ; ttt tann 

fi/2T c (fi/2T c ) 2 2T C 



2 n 
1 H — arctan 

7T TVT, 



(A.9) 



which has the well-known BCS limit: 



(A.10) 
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